BOUNDARY REGULARITY OF ROTATING VORTEX PATCHES 
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Abstract. We show that the boundary of a rotating vortex patch (or V-state, in the 
terminology of Deem and Zabusky) is C°° provided the patch is close to the bifurcation 
circle in the Lipschitz norm. The rotating patch is also convex if it is close to the 
bifurcation circle in the C 2 norm. Our proof is based on Burbea's approach to F-states. 
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O! 1. Introduction 

(N 

The motion of a two dimensional inviscid incompressible fluid is governed by Euler 
equations 

( d t v(z,t) + (v- Vv)(z,t) = -Vp(z,t), zeC, t > 0, 
(1) I divt» = 0, 

v(z,0) = v (z), 

where v(z, t) is the velocity field at the point (z, t) G C x E + and p is the pressure, which 
is a scalar function. The operators v ■ V and div are defined by 

v ■ V = Vidi + t> 2 <9 2 and divt> = diVi + 92^2- 

The velocity field v is divergence free because the fluid is incompressible. We recall that 
in dimension two the vorticity is given by the scalar u> = d\V2 — d^v \ and it is easy to 
see that we can recover the velocity through u by using the Biot-Savart law. Indeed, 
identifying v = (t>i,t>2) with v\ + if 2 and performing a simple calculation one gets 

2dv = iu, with d := d z = -(<9i — id 2 ). 
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Since z i— > -4= is the fundamental solution of the complex operator d we get the Biot- 
Savart law 

(2) v(z,t) = ^ I ^-dA((), zeC, 

2tt J c z- C, 

dA being planar Lebesgue measure. Taking curl in the first equation of the system (CQ) 
one obtains the vorticity equation 



(3) 



d t oj + v ■ Vw = 0, 
u(z, 0) = u (z), 

where uq denotes the initial vorticity and v is given by the Biot-Savart law ([2]). Equation 
([3]) simply means that the vorticity is constant along particle trajectories. Under mild 
smoothness assumptions the Euler system is equivalent to the vorticity formulation Q- 
([2]). A convenient reference for these results is |BM[ Chapter 2]. 

It is a deep fact, known as Yudovich Theorem, that the vorticity equation has a unique 
global solution in the weak sense when the initial vorticity Uq lies in L l D L°°. See for 
instance |BM[ Chapter 8]). A vortex patch is the solution of ([3]) with initial condition 
the characteristic function of a bounded domain Dq. Since the vorticity is transported 
along trajectories, we get that u(z,t) = XD t { z )i where D t = X(D ,t) is the image of D 
by the flow. Recall that the flow X is the solution of the nonlinear integral equation 

X(z,t) = z+[ v{X(z,T),T)dr, zeC,t>0. 



In the special case where D is the open unit disc the vorticity is radial and thus we get 
a steady flow. In particular, D t = D , t > 0, and the particle trajectories are circles 
centered at the origin. A remarkable fact discovered by Kirchoff is that when the initial 
condition is the characteristic function of an ellipse centered at the origin, then the domain 
D t is a rotation of Dq. Indeed, D t = e ltn Do, where the angular velocity Q is determined 
by the semi-axis a and b of the initial ellipse through the formula Q = ab/(a + b) 2 . See 
for instance [BMl P-304] and |Q p.232]. 

A rotating vortex patch or V-state is a domain Dq such that if xd is the initial 
condition of the vorticity equation, then the region of vorticity 1 rotates with constant 
angular velocity around its center of mass, which we assume to be the origin. In other 
words, D t = e ltQ D or, which is the same, the vorticity at time t is given by 

u(z,t) = X D (e- ltn z), zeC, t>0. 

Here the angular velocity Q is a positive number associated to D . 

To the best of our knowledge the ellipse is the only V-state for which a closed formula is 
known. Deem and Zabusky wrote an equation for the V-states and solved it numerically. 
They found V-states with m-fold symmetry for each integer m > 2. It is extremely 
interesting to look at their pictures in [DZJ : one can see domains with smooth boundary, 
which evolve with certain associated parameters to produce in the limit domains whose 
boundaries have corners at right angles (see |WOZ] ). 

Burbea gave a mathematical proof of the existence of m-fold symmetric V-states using 
bifurcation from the circle solution [B]. His approach is, in Aref's words, elegant and deep 
[A~l p. 346]. It is indeed so. He finds an equation for the V-states and then uses conformal 
mapping to rewrite the equation in a functional analytic framework in which bifurcation 
theory can be applied. Unfortunately the proof has a gap, which occurs when the space 
to which bifurcation theory should be applied is set up. The suggestion made in [HI p. 8] 
of using the standard Hardy space does not work. One reason is that the operator Q(f) 



in [Bj p. 8] involves one derivative of / and functions in the Hardy space do not necessarily 
have derivatives. Another reason is that one needs a space which guarantees that small 
analytic perturbations of the identity in the space are conformal. The reader will find 
an exposition of Burbea's approach in section 2 and a complete proof of the existence 
of m- fold symmetric V-states in section 3. It is our impression that this beautiful and 
striking theorem deserves to be more widely known than it appears to be now. In section 
4 we prove our main result, which states that if a bifurcated V-state is close enough to 
the circle in the Lipschitz sense then its boundary is of class C°°. We also show that a 
bifurcated V-state is convex if it is close enough to the circle in the C 2 norm. There is a 
dark zone, where the V-state has boundary of class C l+a for some positive a but it is not 
close enough to the circle, in which the actual smoothness properties of the boundary are 
unknown. Also the nature of the singularities of the "limiting" V-states of Wu, Overman 
and Zabusky ( |WOZj ) are not well understood at all. 

We adhere to the convention of denoting by C a constant independent of the relevant 
parameters under consideration. The constant may change its actual value at different 
occurrences. 



2. Burbea's approach to V-states 

We first derive an equation for simply connected vortex patches which have smooth 
boundary for all times. Consider two parametrizations of the boundary 8D t of the patch 
at time t, say z(a,t) and r]((3,t), and assume that they are proper in the sense that they 
establish a homeomorphism between the interval of definition of the parameters with the 
extremes identified and dD t . Assume also that they are continuously different iable as 
functions of the parameter and time. Then there exists a change of parameters a({3,t) 
such that rj(/3,t) = z(a((3,t),t) for all f3 and t and so we have 

| (ftt) ^ (M) ^) + | (M . 

Since ^-(a,t) is a tangent vector at the boundary at the point z(a,t) and ^((3,t) is a 
scalar we conclude that 

dti Oz 

(4) -^(P,t)-n=—(a,t)-n 

where n is the exterior unit normal vector at the point z(a,t) = rj(/3,t) and the dot 
stands for scalar product in M 2 = C. Thus the quantity in fll]) does not depend on the 
parametrization chosen and represents the speed of the boundary in the normal direction. 
On the other hand v(z(a,t),t) ■ n is the normal component of the velocity of a particle 
which is at the point z(a, t) at time t. Since the boundary is advected by the velocity v 
we get the equation 

Oz 

(5) —(a,t)-n = v(z(a,t),t)-n 

which describes the motion of the boundary of the patch. 
Now we introduce the stream function 







(i logl 


■\*XdA 



(6) il>(z,t) = (^ lo sl ■\*XD t )(z 

Clearly <91og \ z\ 2 = 1/z, where d = d/dz. Hence 

2idif)(z,t) = ^- *XD t ) {-) rl i.t) 



and 

(v ■ n)(z(a,t),t) = (V" 1 ^ • n)(z(a,t),t) 
= -{V^-r){z{a,t),t) 
dip 

= --j^{z{a,t),t). 

where f is the unit tangent vector and s is the arc-length parameter of the curve dD t . 
Therefore ([5]) becomes 

dip dz{a,t) _ 

(7) Ts {z ^ t)) = — dT- n - 

Assume now that the patch rotates with angular velocity Q. Let Zq{oP) be a proper 
continuously differentiable parametrization of 8Dq and set z(a,t) = e tQt zo(a), which is a 
proper parametrization of dD t . Then 

dz 

— (a,t) = iQz(a,t) 

and (j7|) becomes 

(8) — (z(a, t),t) = —iQz(a, t) ■ n = Qz(a, t) ■ f. 

ds 

Taking a = s, the arc-length parameter on dD , we obtain 

■^-\z(s,t)\ 2 = 2Re(z(s,t)T) = 2z(s,t) ■ f 
ds 



so that, by © 



dip . . . . Vt d . . . l9 



and integrating with respect to s yields, for a certain constant c(t) depending on t, 

(9) iP(z,t) = ^\z\ 2 + c(t), zedD t . 

Since the steps can be reversed this is the equation of V-states. 

The goal now is to use conformal mappings to "parametrize" V-states. For that one 
needs to modify the preceding equation to get a form more amenable to the use of 
analytic functions. Fix t and take derivatives in with respect to s on dD t . We just 
get a restatement of (jSJ), that is, 



(10) 2 Re yr^{z{ s i *K 0> = Re (nz{s, t)z'(s, t) 

where the prime means derivative with respect to s. By the generalized Cauchy formula 
(which follows from a direct application of Green-Stokes) one has 

Taking the d = d/dz derivative in ([6]) and applying (fTTjl 



dz n J D 



~~ -2 ~ o — r / t «C> z E Dp 



z-( 2m J dDt C-z 



The first identity above implies that z i-)- ^(z,t) extends continuously to the closed 
domain D t , since the Cauchy integral of a bounded compactly supported function is quasi 
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Lipschitz (its modulus of continuity is 0(8\ \og5\), 5 < 1/2). Thus the same happens to 
the Cauchy integral of the function ( on dD t . Hence (fTUj) becomes, with A = 1 — 2Q, 

Re(\zz'(s,t)- -L ! J^dCz'{s,t))=U, zedD t , 

\ 27T2 J dDt (, — Z J 

where the integral over dD t for z G dD t should be understood as the limit as w G D t 
tends to z of the corresponding integral for w. 

Integrating with respect to s on dD t we conclude that, for some constant c(t) depending 
on t, 

(12) \\z\ 2 + 2Re^- [ Clog f 1 - 7 ) d( = c(t), z G dD t . 

2m JdD t V C/ 

Remember that the origin belongs to D t and thus, for each C G dD t , there exist a branch 
of the logarithm of 1 — z/( in D t \ {(} taking the value at z = 0. Also fixed z G dD t 
there exists a branch of the logarithm of 1 — z/Q in \ D t U {2} taking the value at 
C = 00. Therefore the integral in ffT2l) exists for each z G dD t and defines a continuous 
function on <9_Dt. Later on we will differentiate this integral with respect to z on dD t . 

Notice that the equation ffl2l is invariant by rotations and dilations. Consequently, 
using the fact D t = e ltn D and performing a change of variables equation f[T2]) reduces to 

,2 — 1 /" ^ z 



(13) A|z| 2 + 2Re— / Clog I-7 UC = c, z e dD . 

We have used the fact c(t) does not depend on t since the left hand side of f[T5|) is inde- 
pendent on the time variable. Equation f lTB"]) is Burbea's equation for simply connected 
V-states. 

We introduce now conformal mappings. Let Ee> = \ D be the exterior of D Q and 
Ea = Coo \ A the exterior of the unit disc A = {z : \z\ < 1}. Let $ be a conformal 
mapping of E& onto -Ed preserving the point at 00. This mapping can be expanded in 
Ea as 

$(z) =a(z + Y, 

n>0 

for some complex number a. Making a rotation in the independent variable z we can 
assume that a is a positive number and then dilating the domain we can further assume 
that a = 1. Following Burbea we change variables in ffT3|) setting z = $(u>) and C = 3>(t). 
Then ffTBl becomes 



«7 
Z r ' 



(14) 



A|$(w)| 2 + 2Re— /" $(r) log ( 1 - ] $'(r) rfr = c, wGl 



Recall that we are assuming <9-Do to have rectifiable boundary. It is then a classical result 
that $ can be extended continuously to the unit cercle T = {,zgC:|2:| = 1} and the 
extension is absolutely continuous, so that <3>'(u>) exists for almost all w G T and the 
resulting function is in L X (T) (see [P]). Later on we will work with a $ whose extension 
to T is of class C 1+a (T) so that $' will be Holder continuous of order a on T. To simplify 



the notation, set 



(15) a($)H = ^^$(r)log(l-|^)$'(r)rfr, »GT 



and 



(16) m($,A) = — /(A|$(w)| 2 + 2Recr($)(w)) \dw\, 

2% J j 



where \dw\ denotes the length measure on the circle. Hence (I14p is 

(17) X\$(w)\ 2 + 2Rea($)(w) - m($, A) = 0, w E T. 

It is more convenient to set <&(z) = z + f(z) with / analytic on and define the operator 

(18) S{f)(w) = <r(*)(v)). 

That this is a good point of view is confirmed by the fact that 

(19) S(0)(to) = -^ /log (l --) — = (), wgT, 

because the integrand is analytic on and has a double zero at oo. Thus ()17p is satisfied 
for $(2) = 2 and each A. Now define 

(20) F(X, f)(w) := X\w + f(w)\ 2 + 2 ReS(f)(w) - m(I + /, A), weT, 

where I stands for the identity function. Clearly Burbea's equation can be rewritten as 

F(A,/) = 0. 

We know that 

F(A,0)=0, XeR, 

which simply says that the disc rotates with any angular velocity. Burbea's idea at this 
point is to apply bifurcation theory in order to prove the existence of m-fold V-states. We 
will apply in the next section Crandall-Rabinowitz's Theorem whose original statement 
in \CR\ p. 325] reads as follows. 

Crandall-Rabinowitz's Theorem. Let X , Y be Banach spaces, V a neighborhood of 
in X and 

F: (-1,1) xV 

have the properties 

(a) F(t,0) = for any \t\ < 1. 

(b) The partial derivatives F t , F x and F tx exist and are continuous. 

(c) ^(^(0,0)) and Y/ R(F x (0, 0)) are one- dimensional. 

(d) F tx (0,0)x i R(F x (0,0)), where 

N{F x {0,0)) = span{x }. 

If Z is any complement of N(F x (0, 0)) in X, then there is a neighborhood U of (0,0) in 
Rxl, an interval (—a, a), and continuous functions (p: (—a, a) — > R, ip: (—a, a) — > Z 
such that <p(0) = 0, ip(0) = and 

F-\0) HU = {(</?(£), £z + ■ ICI < a} U {(t, 0) : (t, 0) G U}. 

3. Existence of m-FOLD V-states 

In this section we will apply Crandall-Rabinowitz's Theorem to prove the existence of 
m-fold V-states for each integer m > 2. For m = 2 one recovers the Kirchoff ellipses. 
As a by-product of this formalism we get a Holderian boundary regularity result for the 
V-states which are close to a point of the bifurcation set {(A, /) = (— , 0), m = 2, 3, . . . }. 
Later we will see how to establish the C°° regularity of the boundary of these V-states 
by using hidden smoothing effects of the nonlinear equation (JXTJ) . First, we establish the 
following result. 
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Theorem 1. Given < a < 1 and m = 2,3,... there exists a curve of m-fold rotating 
vortex patches with boundary of class C 1+a bifurcating at the circle solution. 

More precisely, there exist a > and continuous functions A: (—a, a) — > R ? ip: (—a, a) — > 
C 1+a (T) satisfying A(0) = 1/m, ip(0) = 0, such that the Fourier series of is of the 
form 

VK£)H = a^Ow 2 ™- 1 + ■■■ + a nm ^{i)w nm - 1 + ..., WET, 

and 

^Wa^ + ^fOH) = 0, f e (-a, a). 

The mapping 

= z{i + ^ + e« 2m -i(0 ^ + • • • + Mfl + • • • ) 

zs conformal and of class C l+a on C\ A, and the complement o/$^(C\ A) is an m-fold 
rotating vortex patch. 

The proof of this theorem requires some lengthly work, which will be presented in 
several steps. We start by introducing the spaces X and Y and then we will check all the 
assumptions of Crandall-Rabinowitz's Theorem. 

3.1. The spaces X and Y . The choice of the spaces X and Y is a key point, which was 
overlooked in [B] . Before giving the complete description of these spaces we need first to 
recall the definition of the Holder spaces C n+a (Q). Let Q be a nonempty open set of M. d 
and < a < 1. We denote by C a (Q) the space of continuous functions / such that 

i, fl i 1/0*0- /(y)l r 

C-ip.) ■= / L°° + SUp : < 00, 

x ^yen \x - y\ a 

where ||/||l°° stands for the supremum norm of / on Q. More generally, for a non-negative 
integer n the Holder space C n+a (fl) consists of those functions of class C n whose n— th 
order derivatives are Holder continuous with exponent a in Q. It is equipped with the 
norm 

||/||o»+«(n) = ll/IU- + Yl W dl fWc a (n)- 

\y\=n 

We will also make use of the space C l+a (T) which is the set of continuously differentiable 
functions / on the unit circle T whose derivatives satisfy a Holder condition of order a, 
endowed with the norm 

df 

||/||cfi+«»(T) = ll/IU- + ~j- 
where || • || a is the usual Lipschitz semi- norm of order a 

H H \g(x)-g(y)\ 

\\9\\a = sup — ; — 

x^yaf \x-y\ a 

We define in a similar way the spaces C n+a (T), for each positive integer n and a e]0, 1[. 
A word on the operator d/dw is in order. Any function / : T — > K can be identified with 
a 27T— periodic function g : R — > M via the formula 

f(w) = g(6), w = e ie . 

Therefore for a smooth function / we get 

-ie-V(f). 



dw 



It will be more convenient in the sequel to work with d/dw instead of d/d9. Since they 
differ only by a smooth factor it really makes no difference. Notice that we have the 
identity 

(2i) 131 = - _LE. 

dw w 2 dw 

On the other hand, if we denote by C^ a (M.) the subspace of C 1+Q (K) consisting of 2n- 
periodic functions, then we can identify C 1+a (T) with C^^M) and 



Let C\ +a (A c ) be the space of analytic functions on \ A whose derivatives satisfy a 
Holder condition of order a up to T. This is also the space of functions in C 1+a (T) whose 
Fourier coefficients of positive frequency vanish. 
The space X is defined as 

oo 

(22) X = |/ G C 1+a (T);f(w) = ^a n «J", w G T, a n G R, n > C)} 

n=0 

and coincides with the subspace of C^ +a (A c ) consisting of those functions in C* +Q (A C ) 
whose boundary values have real Fourier coefficients. Later on we will modify appropri- 
ately the space X to get m-fold symmetry, but for now it is simpler and clearer to work 
with this X. The reader will understand later why we require the a n to be real. For the 
time being we just comment on the geometric meaning of this requirement. If we set, 
for / G X, &(z) = z + f(z) and $ happens to be conformal on \ A, then the com- 
plement in Cqc of the closure of $(Coo \ A) is a simply connected domain D symmetric 
with respect to the real axis. Conversely, if one starts with a bounded simply connected 
domain D symmetric with respect to the real axis and $ is the conformal mapping of the 
complement of the closed unit disc onto the complement of D, then the coefficients in the 
expansion of $ as a power series in \ jz are real. Now, if one is given a domain D with 
an axis of symmetry containing the origin, after a rotation one can assume that this axis 
is the real line. Domains with m-fold symmetry have an axis of symmetry, so it is not 
really restrictive for our purposes to work with functions with real Fourier coefficients. 

Let V stand for the open ball with center and radius 1 in C^ +a (A c ). If / is in V 
then the function = z + f(z) is analytic on {z : \z\ > 1) and is injective there. For, 
by the maximum principle, 



df 



dw 



Mi]) { : |*| > 1, \w\ >l,z^w} := <) < I. 



L°°(A C ) ' \Z — W\ 



and hence 

- ${w)\ >\z-w\- \f{z) - f(w)\ >{1-S)\z- w\, \z\ > 1, \w\ > 1. 

For / G V define S(f) as in ffTSl) and (TT51) . where $(z) = z + f(z). The definition makes 
sense precisely because $ is injective and thus a branch of the logarithm of 1 — $>(w) /$(t) 
can be defined taking the value 1 at oo (that is, as r — > oo), as we argued after (Tl3i) . We 
define now a function F(X, /) on R x V by 

(23) F(X, f)(w) = X\w + f(w)\ 2 + 2 Re S(f)(w) - m(Id + /, A), 

where S(f) is as in (ITS]) and m(Id + /, A) as in (IT6]) . This is the function F to which we 
will apply Crandall-Rabinowitz's Theorem. 
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We define now the space Y as the subspace of C 1+a (T) consisting of real-valued func- 
tions with zero integral and real Fourier coefficients. More precisely, 

(24) Y := G C 1+a (T) : g(w) = ^ a n w n , w6T and a n = a_ n G R,n > o|. 

Since we have subtracted the mean in (123]) it is clear that F(X,f) is a real- valued 
function with zero integral. To show that F maps X into Y it remains to show that 
F(X, f) belongs to C 1+Q (T) and that its Fourier coefficients are real. 

3.2. F(X, /) is in C 1+Q (T). Since F(X, f) has zero integral its norm in the space C 1+a (T) 
is equivalent to the C a (T) norm of its derivative w t— > j-{X, f)(w). Now, since C 1+a (T) 

is an algebra the problem reduces to showing that w i— >■ (w) belongs to C a (T). For 
that we need to compute the derivative of S(f)(w) with respect to w. 

Lemma 1. For any bilipschitz function $ : T — > C of class C 1 (T) we have 



(25) ^ a)W ._ 4>) ^ + ,,^/ T _|i)^ dT j, raeT . 

Proof. The proof of this lemma consists in computing the derivative with respect to w in 
the sense of distributions by integrating against a test function on T. When integrating 
by parts we will get —<&'(«?) times the principal value integral in (125]) plus a "boundary 
term" — $'(u>)$(u>)/2, which is due to the fact that the logarithm in the definition of 
S(f) is not continuous on the diagonal. 

An alternative argument goes as follows. Bring the integral defining S(f) on the C 1 
Jordan curve T = $(T) to obtain the function 

o-(z) = ^- [ c Mi - j) dc, ^r. 

2m J r Q 

Think of a as an analytic function of z G D, the domain enclosed by T. Its derivative is 

i r 7 l f l 



Z-ni JvC- z ^ JdC 

where the identity is the generalized Cauchy formula for the function z. The integral 
on D in the right-hand side is a quasi-Lipschitz function, as the Cauchy integral of a 
bounded function, and hence the left-hand side extends continuously up to the boundary. 
Call this extension —C(Q(z), z 6 T, where the notation refers to "boundary values of 
the Cauchy integral of the function ( " . By Plemelj's formula (see, for example, (Vj p. 
143]) 

(26) C(C)(*) = |+p.v. ± [J-dC. 

We conclude that a is differentiable on F and its derivative with respect to z is given by 
minus the expression in f l2T)j) . Changing variables to return to the unit circle and applying 
the chain rule we get (|25p . □ 

Now everything is reduced to checking that the the principal value integral in f[2"5"j) 
satisfies a Holder condition of order a. This can be done rather easily in at least two 
ways. The first consists in considering the operator 

(27) Tg(w) = p.v. J- [ — p-L-. - dr, w G T, 

9 



and showing that T maps boundedly C Q (T) into itself. This can be achieved by applying 
the T(l)— Theorem for Holder spaces of Wittmann [Wl Theorem 2.1, p. 584] (see also 
[Gaj ). Before stating Witmann's result recall that the maximal singular integral of g is 



T*(g)(w) 



sup 

e>0 



9{r 



dr 



we J, 



>r\D(w,e) 

D(w,e) being the disc centered at w of radius e. Wittmann's T(l)— Theorem asserts 
in our context that the C a (T) boundedness of T follows by checking that the kernel 
is "standard", that T*(l) is bounded and that the operator T applied to the constant 
function 1 lies in C a (T). That the kernel is standard means in the situation we are 
considering that 



(28) 
and 
(29) 



1 



d 



1 



< 



C 



r — w\ 



-, r, w G T 



dw \ $(r) - $(u0 



< 



C 



r, w G T 



which are clearly satisfied because $ is bilipschitz on its domain. To compute T(l)(w),w G 
T, denote by j e , e > 0, the arc which is the intersection of the circle centered at w of ra- 
dius e and the complement of the open unit disc, with the counter-clockwise orientation. 
Let T e be the closed Jordan curve consisting of the arc 7 e followed by the part of the unit 
circle at distance from w not less than e. We claim that 

, If t — w dr If t — w dr 

T(1)H 



lim 

e-j-0 



(30) 



2ttz 
1 



$(t) - <&(w) 



w 



2-ki 



7e 



T 



w 



2<$>'{w) 



The integral over T e is 1 since the integrand is analytic as a function of r on the exterior 
of the unit disc and we have 



lim 



T — W 



\t\^oo $(t) - $(U>) 

The limit as e — > of the integral over r y e is 

— ; — r lim / = 



1. 



7T? 



and so fl30|) is proven. From the assumption $ G C 1+a (T) combined with f )30|) we obtain 
that T(l) G C a (T). It is also clear that T*(l) is bounded on T from the argument above. 

For future reference we record now the following identity, whose proof is similar to that 
just described of fl30l . 

1 f $'(rl 1 
(31) p.v.— / - 1 ' f dr=-- 

K ' 1 2m Jj $(r) - 2 

The second argument to show that the principal value integral fTSBT) satisfies a Holder 
condition of order a consists in changing variables ( = $(t) and z = $(w) and pass to a 
principal value integral on the C 1+a Jordan curve T = $(T). The kernel of the operator 
one obtains is the Cauchy kernel and the integration is with respect to dQ. This operator 
sends C a (T) into itself on curves satisfying a mild regularity assumption called Ahlfors 
regularity. We omit the details. 
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3.3. Real Fourier coefficients. We intend to show that if / G X then F(\, f) has 
real Fourier coefficients. Let R stand for the space of continuous functions on T with 
real Fourier coefficients. We claim that R is an algebra invariant under conjugation and 
under taking inverses of functions which vanish nowhere on T. Moreover if a function in 
R is continuously differentiable then its derivative stays in R. The only non obvious fact 
is the stability under taking inverses. 

Lemma 2. If f G R and f(e ie ) ^ for all real 6, then 1/f G R. 

Proof. Let f(e %e ) = Y^=-oo a n eme an d (j)( eie ) — Y^=-oo a n^ me be the Fourier series 
expansions of / and j. Set a n = b n + ic n with b n and c n real. Then 



1 oo 



2^ ((a*b)(n)+i(a*c)(n)} 



And 



and thus 

(32) (a * c)(n) = 0, for each integer n. 



oo 

— oo 



Here we are using the notation a * b for the convolution of the sequences a = (a n )" 
and b = (&n)^_oo; tliat is > ( a * b )( n ) = Ylm=-oo a n~mb m - Let 9 i-)> ^(e^) be the func- 
tion with Fourier series X^-oo c ri e ™ 9 - Then g G i^ 2 (T), because the sequence (an)™.^ 
is square summable and so (On^-cc is. By f l3"2"j) f(e ld )g(e l6 ) = for almost all 6. 
Hence g{e ld ) = for almost all 6, which means that c n = 0, for all n. Therefore a n = b n 
is real for all n and 1/ / G R. □ 



From the discussion above it is clear that the term A|$(w)| 2 = X[w + f(w))[w + f(w)] 
has real Fourier coefficients if / does. Thus the following lemma completes the proof that 
F(X, f) has real Fourier coefficients if / G X. 

Lemma 3. If f £ X then ReS(f) has real Fourier coefficients. 

Proof. We are going to show that, in fact, S(f) has real Fourier coefficients. Set 

oo 

71=0 

By Schwarz inequality 

oo oo oo 

(33) J2\ a n\ < (X^ 7 * 2 ! -! 2 ) 172 (5Z^2") 1/2 = C||d//^|U 2(T) < C||d//^|Uoc (T) < OO. 

n=l n=l n=l 

Arguing as we did for (}T9]) we conclude that 

and, consequently, the function S(f) takes the form 



S(f)(w) = ^~l f(r) log (l - S$ ) $'(r) dr. 
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Thus (I3"3"j) yields 



(34) S(f) (w) = > >„-^ / t" log I 1 - |£4 I *'(t) <ir, «, 6 T. 



n=0 

Choose A > 1 big enough so that 

(35) |$(r)| > 2||$|| L o 0(T) , \t\ = A. 

Since the integrand in (13"1"|) is analytic in |r| > 1, the integral there can be taken on 
|t| = A Then we can use the expansion 

to obtain 

OO OO , 

(36) s(/)H = -£«»£n^*(«0 



where 



k 

n=0 A:=l 



A n k = — f r n - fe $'(r) ( -^—] dr 



(37) = i.f ^ fc+1$, ^(^y) fc '^ 



is the Fourier coefficient corresponding to the frequency — (n — k + 1) of the function 
t i— > $'(r)(r/$(r)) fe . Clearly r i— )■ $(r)/r has real Fourier coefficients and thus so does 
r i — y $ / (r)(r/$(r)) fe , by Lemma [2j Hence A n ^ is a real number and so (1361) says that 
S^/) has real Fourier coefficients □ 



3.4. F is Gateaux different iable. We show here that F is Gateaux differentiable and in 
the next subsection we will check that the directional derivatives are continuous. This will 
show that F is continuously differentiable on its domain R x V. In proving differentiability 
properties of F(X, f) we can consider only the first two terms in (T2U|) because the third, 
m(I + /, A), is a mean of the other two. Call G(X, f) the sum of the first two terms in 
(J2QD. 

There is no problem with the partial derivative with respect to A because the depen- 
dence on A is linear. We get 

D x F(X,f)(w) = \w + f(w)\ 2 

which is obviously continuous with respect to / in the topology of X. 

Take now h £ X and compute the derivative of G(X, /) with respect to / in the direction 

h, that is, 

D f G(XJ)(h) = SLG(\,f + th)> 

at 1 1=0 

(38) = 2A Re(<M) + 2 Re ^-S(f + th) i . 

at | t=0 

12 



The computation of ^S(f + th)\ t=0 yields the four terms below, 

+ + ^L W) ' og ( i -W)) h ' {T)dT 



1 f -r-rx h(w) - h(r) 
2tu J t v ; $ (r) - $H 



+ — [ ^QhM&MdT 

(39) := /»)(«;) + /»)(«;) + C7(/, /»)(«;) + !>(/,/»), 

where the last identity is the definition of the functions A, 5, C and .D. Notice that D is 
independent of w. 

We proceed now to prove that A(f,h) G C 1+Q (T). For the estimate of the absolute 
value of A(f, h)(w) we can assume without loss of generality that w — 1. Computing the 
derivative of log(l — $(l)/$(r)) with respect to r in |r| > 1 and using the fundamental 
theorem of calculus we se that 

40 log =-/ — / ',\', r rdt, l^r, r > 1. 

1 ^ V $(^)/ A (®(tr) - $(l))$(tr) II — 

The preceding identity is very useful in estimating the integrals containing a logarithmic 
term, as we will see below. We proceed now to estimate from below the factors in the 
denominator of the fraction inside the integral in (T4*0|) . The function $(r)/r is analytic 
on C \ A and takes the value 1 at oo. Hence, by the maximum principle and recalling 
that / G V, we obtain 

(41) 1 - II/IU-ct) < < 1 + \\fh°°(T), \r\ > 1, 
and 

(42) Mr) - $(1)| > (1 - ||/'|U~( T) )|t - 1|, |r| > 1. 
Therefore, by flU}, flU} and (T42]) we get for r G T\{1}, 



log 1 



$(r) 



< 1^(1)111^11^00^0) 



dt 



< 



\<$>(tr) - $(l)||$(tr)| 

(l + ll/IU-(T))(l+||f||L-(T)) f X dt 



(i - II/IU~(to) (i - Il/'lU-cr)) A t\tT-l\' 

We split the interval of integration in the last integral above in three subintervals : 

(1, 1 + \t — 1|), (1 + |r — 1|, 3) and (3, oo). In the first we notice that for t > 1 and |r| = 1 

we have the elementary inequality 

\tr - 1| = \t - t\ = \t - t\ > |1 - t\ 

and so 

" 1+ |r-H dt 

-, r < !■ 

t\tr - 1| - 

The integral on the second interval can be estimated straightforwardly as follows 

" 3 dt r 3 dt , 2 



-]r-l| ~~ 1| il+|r-l| * 
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T 



Finally 

dt f°° dt ,3 

< / 777 77 = lo S ■ 



t\tr - 1| - J 3 t(t-l) 2' 
Therefore, collecting the preceding inequalities, 

(43) ■ - 



log 1 



$ r 



<C7(/)(l + |log|r-l||), rGT, 



where C(/) is a constant depending only on the C 1 (T) norm of /. 
Integration in r on T readily yields 

The next step is to estimate the uniform norm and the Lipschitz semi-norm of order a 
of the function w i— > -4-A(f, h){w). As in Lemma [T] one has 



(44) -^(/, ft) W = (-U + p . v . _ ^ >( U_ ^ dT | , roe T, 

and the only difficulty lies in estimating the C a (T) norm of the principal value integral. 
But this has already been done in subsection 3.2. Therefore A(f, h) G C 1+a (T). 

The proof that the term B(f, h) belongs to C 1+a (T) is basically the same. The expres- 
sion for the derivative is 



— B(f, h)(w) = -&(w) ' ' + p.v. / , / —, - 

dw V ; I 2&(w) 2m J r $(r) - <$>{w) 

and one deals with the principal value integral as before, using the operator in (J2j 

The proof that the term C(f, h) is in C 1+a (T) contains a little variation of the preceding 
argument. First, we observe that the quotient 

h(r) — h(w) 
$(t) - 

takes continuously the value h' (w) / (w) on the diagonal of T. Consequently, in com- 
puting the derivative of C(f, h){w) no boundary terms will arise and we get 

4-C(f,h)(w) = h'(w)p.v.— [ — 9 ^\ / . dr 

2m y w=1 $(r) - $0) 

+ &(w) p.v. / — ■) / — ^ i— g(T)dr 



where g(r) = $(r) $'(t). The mapping properties of the operator T in f |27|) take care of 
the principal value integral in the first term. We view the principal value integral in the 
second term as an operator U acting on g. Its kernel is standard, because $ is bilipschitz, 
and the action of U on the constant function 1 is, as in ( 130]) . 

U(l)(w) = lim/ . ff? ~ dr 

1 f h(r) — h(w) ( r — w \ dr 
+ hm ' ' 1 



£->0 2m' T — w \ $(t) — / T — w 

h'{w) 
2$W' 
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In the integral on T e the integrand is analytic on C \ A and has a double zero at oo, and 
so the integral vanishes. The limit as e tends to of the integral on r y e is m times the limit 
as r tends to w of the quotients inside, that is, h' \w) / $' (w) 2 . Then U(l) G C a (T). The 
previous discussion also gives that U*(l) is bounded. Thus one can apply Wittmann's 
T(l)— Theorem (see the paragraph after (127)1 ) and conclude that U maps C a (T) into 
itself. This completes the proof that w DfF(X, f)(h)(w) G C 1+a (T). On the other 
hand, it is clear that DfF(X, f)(h) depends linearly on h. Therefore F(X,f) is Gateaux 
differentiate at any point (A, /) G R x V. 



3.5. Continuity of DfF(X, f). In this subsection we prove that the mapping DfF(X, f) 
is continuous as a function of / G V taking values in the space of bounded linear operators 
from X into Y. In particular this shows that F(X, f) is continuously differentiate in the 
Frechet sense. What one has to do is the following. Fix / G V and show an inequality of 
the type 

\\D f F(X,f)(h)-D g F(X,g)(h)\\ci + <* m < C 1+a (/)||/ - ^|| c i + « (T) ||/i|| c i +a(T) 

for g G X close enough to /. Here we denote by C\ +a (f) a constant depending on the 
norm By (I3"8|) and ( 150|) this amounts to prove similar inequalities for A(f, h), B(f, h) 

and C(f, h) in place of DfF(X, /). For instance, 

\\Mfi h ) ~ A (9,h)\\c^(j) < Ci +a (/)||/-p|| c a+«(T)||/l||cn+«(T) 
for g G X close to /. We start with the estimate of the uniform norm of the difference 

(45) ||A(/,/ i )-^,/ i )|| 00 <C 1 (/)||/-^|| cl(T) ||/ i |U 00(T) . 

To prove the uniform estimate above it is enough to consider the point w — 1 in T. Take 
g G V close to / and set ^(r) = r + g(r), |r| > 1. We can easily check that 

A(f,h)(l)-A(g,h)(l) = ^JW)^g{l-^)^'(r)-^(r))dr 

(46) := Ti{f,g){h)+I*tf,g)(h). 
To estimate the first term Xi(f,g)(h) we use 



\Mf,g)(h)\ < 7r\\h\U\f - sf\ 



log 1 



\dr\ 



(47) < C(f)\\h\U\f' ~ 9'Woo. 

To treat the second term X 2 (/, g){h), we use the identity (HOI) which yields 

$(l)$'(tr) $(l)$'(tr) 



/C(r) 



! U$(iT)-$(l))$(fr) (*(frr) - *(l))*(tr) 



rdt, rG T\{1}, 



where 



By elementary algebraic computations one gets 

V ; J x ($(tr) - $(l))$(tr) V ; A ($(tr) - $(l))$(tr) 

+ ^ 1 / (tr 7— . V ^7ttt^7 , T , N rrft 
V Vi (*(^) -$(l))$(tr)^(tr) 

/•°° j^frr) / 1 1 \ 

+ u 7i "^yU^r) -$(i) ~ v(jpj -y(i)) T 1 
■= 

i=i 

Coming back to the identity fj40|) . the first term /Ci(r) takes the form 

MT) = -' (1 >T* (1) l0gfl 



$(1) to V $(r). 
It follows according to ( 113]) that 

|/Ci(r)| < - <7|U~ (T) (l + | log \t - 

As in the proof of (|43j) . we obtain 

|*,(r)| < l»(l)lll«'-«'ll^/" wtr) _ 1 ip(tr)| <fr 

poo y 

< ||/' - <7'|U~(T) (l -H I log |x - 1| l) . 

Concerning the third term, we write 

r°° i 

|K.W| < I^DIII*'!!^,!!*-*!!^,/ mr) _ m{tT)mtT / r 

< C 1 (/)||/-g|| L « (T )(l + |log|T-l||). 



To treat the last term /C4 we use 

- $}(tr) -{*- $}(!) 



$(tr)-$(l) *(tr)-*(l) ($(tr) - $(l))(tf(tr) - 

< 4 JI^-^IU»(A«)|*r-l| 



< 4 



|tr-ip 
II/' - 9'\\l°°(Y) 



\tr-l\ 

Consequently, 

Mr) I < Ci(/)||/' - </|U~ (T) (l + I lQ g k " !| 
Therefore, 

|/C(r)| < C 1 (/)||/-^|| C1(T) (1+ |log|r- 1||), r G T\{1}. 
Hence, coming back to the definition of X 2 (/, g)(h) in ( l46p and integrating in r 

IW,»)(/»)l<Ci(/)||/-»||oi m . 
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Finally, 

|^(/,/ i )(i)-^,/ i )(i)|<c 1 (/)||/ i |U c (T) ||/- 5 || cl(x) 

which is ()45p . 

Our next task is to estimate the difference 

A_ mh)iw) -jL A(gMw) 

in C a (T) for g close to / in C 1+a (T). The difficult term in the formula for the derivative 
of the function A(f, h)(w) in (1441) is the principal value integral 



y 27iiJ T $(r) 

i r Hj) -h(w) _,. . , i— - 

2m 7 T $ r - $(w) K ' 2 y ' 



(48) := J(f,h)(w) + -h(w). 

Thus the difference J(/, h)(w) — I(g, h)(w) is 



J(f,h)(w)-J(g,h)(w) = -L / ff r j ^W (r)-fr'(T))rfT 



2m J T $(r) - $(w) 

+ — f (Mr) -Mwf)V'(T)( , : 1 ; : : — : — 7 — flV 

2iriJj yyJ y " V; V$(r)-$(w) *(r)-*(w)/ 

:= ri(/,0)(to)+7S(/, *)(«;). 

Now we need an estimate for 7i (/,<?) and T^i f,g) in C°(T). Both terms have the form 
(49) T X M = / K{w, r)x(r) dr, X e L°°(T) 

where the kernel K(w,t) and the function \ are 



in Ti(f,g)(w) and 



^ w > r J $( r ) -$(«;) *(r) -*(«;)' Xl ' { ) 

in T 2 (f,g)(w). 

Lemma 4. Assume that the kernel of the operator T in (j49p satisfies 

(1) If zs measurable on T x T and 

|%t)| < C , w,rGl 

(2) For eac/i r G T w h-> r) is differentiable in T\{w} and 

C 

\d w K(w,r)\ < | r , w,t6T, w ^ r. 

\W — T\ 

Then 

(50) |TxK)-TxWI<C||x|UooCo(l + lo gl -) \w 1 -w 2 \, Wl ^w 2 eT, 

V \w 1 -w 2 \J 

for some constant C. In particular, T\ G C a (T) for every < a < 1. 
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The proof of this lemma is simple and standard. For the details of a similar result the 
reader is referred to [MOV} p. 419]. 

To get the desired estimates for T\ (/,<?) and 72(/, <?) we just need to check that their 
kernels satisfy the hypothesis of Lemma HI We deal first with Ti(f,g). From (j4"2j) with 
the point 1 G T replaced by an arbitrary toeT, one readily gets 

\K(w,t)\ < C\\h'\\ L «,(p), \d w K(w,r)\ < C\\h'\\ L ^(T)\r - iu| -1 . 

Therefore 

||7i(/,p)||c q (t) < C \\h \\l°°(t)\\^ — $ IU°°(T) 

< c \\h \\L°°(X)\\f ~ gWc 1 ^)- 

To estimate the kernel of %(f,g) we use 
\K(w, t) J 

< C\\h \\l°°(T)\\® ~ ^ IU°°(t)- 
The derivative of K(w,t) with respect to w can also be estimated easily 

\d w K(w,T)\ < \h'(w)' ' ' 



\h(r) - 


-h(w)\ 


($-*)(w)-($-^)(t) 




$(r)- 


${w)\ 


*(t) - *(iw)| 



+ IMt)-aMI 



&(w) V'(w) 



($(r) - $(tu)) (*(r)-*(tw)) 

< C||/i ||x,oo(T)||$ — \& Hl 00 ^) ( 1 + ||$ + ^ ||l°°(T) 



r — w\ 



< Cl(/)II^I|C1(T)II/-^I|C1(T 



r — w\ 



This gives according to Lemma S] 

\\T2(f,g)\\c a (j) < Ci(/)ll^l|ci(T)ll/ — ^llci(T)- 

Hence we get 

IIZ-^II 

and, finally, gathering all previous estimates 

\\Mf: h ) ~ A (9,h)\\ cl+a{T) < C^flWf - g\\ci+a(T)\\h\\ C i+ a (jy 

This concludes the proof of the continuity of the term A(f, h)(w) with respect to /. 

The proof for the terms B(f,h), C(f,h) and D(f,h) given by (1591 follows a similar 
pattern. We omit the details. 

3.6. Kernel and range of DfF(X,0). In this subsection we study the kernel and the 
range of DfF(X, 0). We also find the "eigenvalues", that is, the values of A for which the 
kernel of DfF(X, 0) is not trivial. In fact, the dimension of the kernel for these particular 
values of A turns out to be 1. Then we check that the range has codimension 1 so that 
Crandall-Rabinovitz's Theorem can be applied. 
Letting / = in p@ and (El we obtain 



D f F{X,0)(h)(w) = 2A Re(h{w)w) + 2 Re [ h(r) log(l - -)dr) - m{h), 

V27T2 Jj T / 

where m(h) is the mean on T with respect to \dw\ of the sum of the first two terms in 
the right-hand side above. Hence w i— > DfF(X,0)(h)(w) has zero integral with respect 
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to \dw\ on T. We would like to compute the Fourier series of w i— > DfF(X,0)(h)(w) in 
terms of the Fourier series of h 



CO 



(51) h(w) = ^bnllf 1 , weT. 

n=0 

Since h G X the Fourier coefficients h n of h are real. Using the expansion 

oo 1 

l g(l - -) = - V -W n T n 
n=l 

and computing we get 

co _ co 1 

D f F(X,0)(h)(w) = ]T(A - -)b n ^w n + ]T(A - -)b n ^w n . 



n * — ' n 

n=l n=l 



From the above expression we immediately conclude that the kernel of DfF(X, 0) is non 
trivial only if A = 1/m for some positive integer m. If this is the case, then the kernel is 
one dimensional and is generated by the function w h- > w" 1 ^ 1 = l/w" 1 ^ 1 . It is precisely at 
this point when we use the fact that the Fourier coefficients of the functions in our space 
X are real. Were not this the case we would get a kernel of real dimension 3. Let us now 
look at the range of DfF(X,0) under the assumption that A = 1/m. Clearly 



CO 

(52) D f F(l/m, 0)(h)(w) = 2 > ( ) b n ^ cos(n9) 



w = e w . 



m n 

n=l 



Notice that DfF(l/m, 0)(h) is a function with zero integral and real Fourier coefficients. 
We have shown in subsection 3.4 that it is in C 1+a (T) and thus in the space Y given 
by (j24p . The only Fourier frequency missing in the expansion (|52p is m, so that it looks 
plausible that a complement of the range of DfF(l/m, 0) is the one dimensional subspace 
generated by cos(m#). To prove this we need to show that each g EY with an expansion 
of the form 

CO CO 

g(w)= E 

with real /3 n is equal to DfF(l/m,0)(h) for some h G X. If h is as in (15 ip with real 
Fourier coefficients then the equation DfF(l/m,0)(h) = g is equivalent to 

2(- --)&„-! =P n , 71=1,2,.. 



or, solving for b n 



m n + 1 

bn — — — Pn+l 
2 n + 1 — m 

m m 2 1 

2 2 n + 1 — m 



The solution to D f F(l/m,0)(h) = g is 



2 

/i(uj) = — 7y ) + — w H(w) 



where 

CO 



n=l 
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and 



' n — m 



H(w)= 2^ r 12 ^^ weT - 

n=l, n^m 

The function G is in C 1+a (T), and then in X, because the Cauchy projection 



c n w n H> ^ c n w n 



n=0 



preserves the space C 1+a (T). This is false for the space C 1 (T) of continuously differen- 
tiable functions on T (because the Cauchy projection does not preserve L°°(T)) and this 
is why we cannot choose the space C 1 (T) in the definition of X and Y. It still remains 
to show that H G C 1+a (T), but this is easy. Set 



K(w)= Y — — , weT, 
^ , n-m 

n=l, n-f-m 

so that K G L 2 (T) C L\T) and H = G * K G C 1+a (T). 

To apply Crandall-Rabinowitz's Theorem we still have to check the transversality con- 
dition, that is, that the the second order partial derivative Df\F(l/m,0) of F with 
respect to / and A applied to the function w i— >■ w" 1-1 is not in the range of DfF(l/m, 0). 
Now Df\F(l/m,Q) can be identified with a bounded linear mapping from X into Y. It 
is easy to see that 

D fx F(l/m, 0)(h)(w) = 2Re(h(w)w), h E X. 

Hence 

D fX F(l/m,0)(w m - 1 )(w) = 2Re(w m ) = 2cos(m#), w = e 10 , 

which is not in the range of DfF(l/m, 0). 

Finally one checks without pain that Df X F(\, f)(h)(w) = 2Re(<&(w)h(w)) is a contin- 
uous function on R x V. 

3.7. m-fold symmetry. We showed in the previous subsections how to apply Crandall- 
Rabinowitz's Theorem to the spaces X and Y. The conclusion is that, given a positive 
integer m, we have a continuous curve (Ag, /t) G R x V, defined for £ in some interval 
of the form (1/m — 8, 1/m + 8), such that F(\%, /g) = 0, £ G (1/m — 8, 1/m + 8). Then 
$£(z) — z + fz(z), \z\ > 1, is a conformal mapping of \ A into some domain U% and 
D% = Coo \ U% is a simply connected vortex patch which rotates with angular velocity 
Q$ — (1 — \)/2. We know that -D^ is a domain with boundary of class C 1+a , but nothing 
else can be said about its symmetry properties without further arguments . The m-fold 
symmetry follows by adding a condition to the spaces X and Y. 

Given m, define X m as the subspace of X consisting of those functions / G X with a 
Fourier expansion of the type 

oo 

f(w) = ^a nm . 1 w nm - 1 , weT. 

n=l 

If / is in the open unit ball of X m the expansion of the associated conformal mapping $ 
in {z : \z\ > 1} is given by 

oo 



. z* 

n=l 
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This will provide the m-fold symmetry of the associated patch, via the relation 

<5>{e l27r/m z) = e i27r/m <S>{z), \z\ > I. 

The space Y m is the subspace of Y consisting of those g G Y whose Fourier coefficients 
vanish at frequencies which are not non-zero multiples of m. In other words, the Fourier 
expansion of g is of the type 

oo 

9{ w ) = ^ nm 2 cos ( nm6 0) w = e%(> i 

n=0 

with (3o = and real (5 nm . Notice that, since the generator w w" 1 " 1 of the kernel of 
DfF(l/m, 0) is in X m for m > 2, we still have that the dimension of the kernel is 1. In 
the same way the codimension of the range of DfF(l/m, 0) in Y m is 1. 

However to apply Crandall-Rabinobitz's Theorem to X m and Y m one has to check that 
F(X, /) G Y m if / G X m . This follows rather easily from work we have already done. One 
has to observe that the space A m of continuous functions on T whose Fourier coefficients 
vanish at frequencies which are not integer multiples of m is an algebra, closed in the 
space of continuous functions on T. Next we remark that w i— > <&'(w) G A m . We also 
need the fact that w (-)■ G A m . To show this, set 

g(w) = ^-l, weT, 
w 

so that g G A m and \\g\\oo < 1- Thus 

oo 

An easy computation gives that w i-> A|$(w)| 2 belongs to A m if / G X m . It remains 
to show that w i-> S(f)(w) is in A m . Recall the identities (|36|) and (137|) . First, a n 7^ 
only for indexes of the form n = mq — 1 for some positive integer q. On the other hand, 
A n k is the Fourier coefficient corresponding to the frequency n — k + 1 of the function 
$'(r)(r/$(r)) fe , which is in A m . Hence A n k is non-zero only if n — k -\- 1 — tw for some 
integer r. Therefore the sum in k is only over indexes which are multiples of m. It remains 
to examine at the Fourier coefficients of §(w) k . Now = wg(w) with g G v4 m and so 
$(u>) fc = w k g(w) k is also in A m because only indexes k which are multiples of m have to 
be taken into account. 

Therefore we can apply Crandall-Rabinobitz's Theorem to X m and Y m and finally 
obtain the existence of m-fold symmetric V-states for each integer m > 2. 

3.8. Kirchhoff 's ellipses. For m = 2 we obtain the ellipses parametrized by w G T y 
w + C,w. The real number £ satisfies —1 < £ < 1 and is a parameter which determines 
the shape of the ellipse. The ellipse is centered at 0, has horizontal semi-axis 1 + £, 
and vertical semi-axis 1 — £. The function z 1— >■ z + - is the conformal mapping of the 
exterior on the unit disc onto the exterior of the ellipse. It is instructive to use Crandall- 
Rabinobitz's Theorem to prove that these ellipses rotate. We are going to apply the 
Theorem to the one dimensional spaces X, which is generated by w, and Y, generated by 
w 2 + w 2 = 2 cos(2#), w = e ld . Notice that X is the kernel of DfF(l/2, 0), which then has 
range {0} of codimension 1 in Y. Of course we have to check that F(X, f) sends X into 
Y. Take f(w) = £wJ with |£| < 1, so that <&(«;) = w + £uT. The term 

\§(w)\ 2 = l + £ 2 + tw 2 + £w 2 

21 



is correct because the constant 1 + £ 2 will disappear when subtracting the mean. We can 
compute explicitly S(f)(w) using ( 13 6 p and the remark that the sum in n and k may be 
reversed because the sum in n is finite. We obtain 

The only term that survives is that corresponding to the index k = 2 and the result of 
the integral is £. Thus 

S(f)(w) = -^(w + Zw) 2 

and 

2Re S(f)(w) = - (^(l + e)(w 2 + w 2 )) + 2e 

Again the constant term will disappear when subtracting the mean and we conclude 
that F(X, f)eY. 

A final remark is that, strictly speaking, the conclusion of Crandall-Rabinowitz's The- 
orem is that for some little interval of £ centered at the associated ellipse rotates. But, 
of course, that any ellipse satisfies Burbea's equation ( TT3]) can be proved directly. It is 
interesting to notice that, in this example, the dependence on both z and £ provided by 
the Theorem is analytic. We do not know how general this fact is. 



4. Boundary smoothness of rotating vortex patches 

In this section we prove our main result, namely that if the bifurcated patch is close 
enough to the circle where bifurcation takes place, then the boundary of the patch is of 
class C°°. Before stating the result more formally we remind the reader of the big picture. 
We called V the set of functions in the unit ball of C 1+Q (T) with real Fourier coefficients 
living only at negative frequencies. Each / e V determines a conformal mapping §(z) = 
z + f(z) of the complement of the open unit disc A into some domain containing the 
point at oo. The boundary of the simply connected domain D = C \ $(C \ A) is the 
Jordan curve $(T) and so, since $ G C 1+a (T), the boundary of D is a Jordan curve of 
class C 1+a . Burbea's existence Theorem asserts that for each integer m > 2 there exists 
a small positive number a and a continuous curve /(£), —a < £ < a, taking values in 
V such that the simply connected domain associated with /(£) is an m-fold rotating 
vortex patch. Since /(0) = 0, Do is the open unit disc and one should think that is a 
domain close to the disc, for small values of £ 6 (—a, a), in the topology determined by 
C 1+a (T) . We claim that if is close enough to the disc in the topology given by C 1 (T) 
then the boundary of is of class C°° . Later on we will show that if is close enough 
to the disc in the topology given by C 2 (T) then is also convex. 

Theorem 2. For each integer m > 3, there exists a small positive €q = eo(m) such that 
if f 6 V defines an m-fold V-state D and WfWc^tr) < e o> then D has boundary of class 
C°°. 

Outline of the proof. Before plunging into the details we present a sketch of the proof of 
the Theorem |2j Burbea's equation for V^-states is F(X, f)(w) = 0, w6T, with 

F(X, f)(w) = X\w + f(w)\ 2 + 2 ReS(f)(w) - m(Id + /, A), w e T, 

where S(f) is given by ([18]) and m(Id + /, A) by (TT6]) . Recall that the reason to subtract 
m(Id + /, A) is that the integral of F(X, f) over T be zero. 
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In the previous section we have used bifurcation theory to prove that, given an integer 
m > 2 and < a < 1, there exists a curve of ^-states passing through (1/m, 0) and taking 
values in a little neighborhood of (1/m, 0) in (0, oo) x V. Since V is contained in C l+a (T) 
we conclude that the non-trivial ^-states we have found have boundary of class C 1+a (T). 
The same approach can be adapted with slight modifications to the space C n+a (T), for 
each positive integer n and each < a < 1. This provides curves of solutions with 
boundaries of class C n+a . However the neighborhood of (1/m, 0) containing the curve 
of solutions decreases as n increases and so the C°° the regularity of boundary cannot 
be reached by using Crandall-Rabinovitz as a black box. Since Crandall-Rabinovitz 
depends essentially on the implicit function theorem, we thought of resorting to Nash- 
Moser implicit function theorem for C°°(T). Unfortunately we were not able to implement 
this idea. We realized later that a simpler method works. The idea is to differentiate the 
equation F(X, f)(w) = with respect to w and carefully study the resulting equation. 
We find a surprising smoothing effect for the unit tangent field to the curve $(T) which 
induces in turn a global smoothing effect for the conformal mapping $. To be more 
precise, we compute dF ^^ , which yields a formula for the quotient 



of the form 



(53) q(w) = w , wel 



A)$(tu) + h(w) 



where l\ is the integral 



(54) h(w) = ± [ t^-t (W U'(r)dr. 



2ni J T $(r) - <&(w) 

As we know, a priori $ G C 1+a (T) and thus q is only in C a (T). But fl53l) suggests that q 
might be of class C 1+a (T) provided h(w) is. It is not difficult to compute the derivative of 
h and check that it is in C P {T), < < a. Thus we get that q G C 1+/3 (T), < < a, 
if the denominator in flo^l) does not vanish in T. This is guaranteed by the smallness 
condition ||/||cinr) < e o- Now the smoothness of q is the same as that of the unit tangent 
vector to the curve $(T) and classical results on the smoothness of conformal mappings 
yield that $ G C 2+/3 (T), < f3 < a. One can then iterate the argument and show that 
under the assumption $ G C 2+/3 (T), < f3 < a the function I\ can be differentiated 
twice with respect to w and the second derivative is in C^(T), < (3 < a. Thus $ G 
C 3+l3 (T), < (3 < a. Since the iteration can be performed any number of times we 
conclude that $ G C°°(T). We begin now with the details of the proof. □ 

Proof of Theorem^ To get (|53i) we differentiate the equation F(X,f)(w) = with re- 
spect to w and use (|2T|) . We get 



\(&(w)$(w) - ^'(w)^(w)) + - l.^Ml(w) = 0. 

V w 2 / dw w 2 dw 

Recall that the parameter A is taken in the interval ]0, 1[. According to fl25|) and fl3B . 

:= -$'(w)^(mj) + /i(w)), to6l 



Putting together the two preceding identities and setting q(w) := §t|^ we obtain ( l53|) . 
Let us show that the denominator in (1331) does not vanish on T if H/Hc^t) is small enough. 
Since = w + f(w), w6T, the denominator in (|53|) is 



(55) £>(w) = (1 - A) (w + f(w)) + h{w), w G T. 

Now 



v ; 2th 7 T $(r) - v ; 2ttz 7 t $(r) - v ; 

11 f r — w T , , dr 



7 T $(r) - r 



27TZ 

where Ji is a notation for the second term and we have used the identity 

T — W 

(56) r — w = , r, w G T. 

TW 

On the other hand, by Lebesgue dominated convergence Theorem we have 



2-ki J T <3>(t) - $(u>) r s^o 2iri J\ T \ =1+£ <&(r) - $(w) r 

The integral on the circle of radius 1 + e is 1 for each positive e, because the integrand is 
an analytic function of r in the exterior of the unit disc whose first term in the expansion 
at oo is 1/t. Hence 

h(w) = — w + Ji(w), w G T. 
Plugging this identity into ( 155]) 



D{w) = -Aw + (1 - A)/H + JiO), «J6l 
The estimate of the L°° norm of Ji can be easily performed as follows 



Therefore 



Kl||i«CT) < — |$(r)-$( W )| H $ 

im T ^«, | T _ W | 

< H^'llym (i + ||/-|u„). 

1 - ||/ ||l°°(t) 



1-0(^)1 > A - (1 - A)||/||ioo (T) - ||Ji||i«»(T) 

l + ll/'lloo 



> A-II/IU-II/' 



OO - 



> A 



> k 



l-ll/'l 

C 1 (T) 



where the last inequality holds provided ||/||ci(T) < 4^x- 

Let us now prove that the function idgTh -^i(w) is more regular than one expects 
and belongs to the space G C 1+(3 (T), for any < (5 < a. Since the quotient w §fe\E§fe] 
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extends continuously to the diagonal of T, in taking derivatives inside the integral defining 
h{w) no "boundary terms" will appear. Then it follows from (1211) and (1311) that 

—±(w) = —z&(w) p.v. / ——±-L—dT 

dw y ' w 2 y ' ^ 2ni J T $(r) - 

(57) + $'M P-v. — I 3 T } ~^ W L &(r)dT 
1 ^ 1 ' 2m J T ($(r) - $H) 2 1 ; 

(58) = + nw) p . v . ± [ *V-*W &{T)dT. 
K ' 2w 2 K J y 2m Jj ($(r) - $(w)) 2 V ; 

To make the appropriate Holder estimate we are looking for it is convenient to write 
the principal value integral above as the sum of two terms by adding and subtracting 
<&'(w)(t — w) in the numerator of the fraction. We get 



P - V - 2m A ($(r) - lrjdr 2ttz J t ($(r) - *(u;)) 2 * [T)dT 



(59) + $'(w) p.v. / — — — — &( T )dT 

1 ; K >y 2mJ T ($(r) - $H) 2 1 ; 

Call I 2 (w) the first term in the right hand side above, namely, 



T — W 



2m J T ($(r) - $(w)) 2 

Let us compute the principal value integral in the second term of (1591) . By (I56p 

(60) 

If t-w * / \ i 1 1/"/ t-w \ 2 $'(r) dr 

P- v - o '■ / 7^73 ITt NN9 $ ( r )» r = P-v- 



27ri y T ($(r) - $(w)) 2 to 27rz y T V$(r) - <&{w)J r r - u> 

We compute the principal value integral above by the method used in dealing with ( 130]) . 
Denote by j e , e > 0, the arc which is the intersection of the circle centered at w of radius 
e and the complement of the open unit disc, with the counter-clockwise orientation. Let 
T e the closed Jordan curve consisting of the arc 7 e followed by the part of the unit circle 
at distance from w not less than e, traversed counterclockwise. Then the principal value 
in (EOD is 

2 $'(t) dr 1 r / t-w \ 2 &(t) dr 

2~i .1 ■ \ <I>( t) <!>( f/'i /' t t-w e-x> 27TZ 7 Te V$(r) - r r - w 

1 /" / t-w \ 2 $'(r) rfr 
2ni J V$(r) - r t-w' 

The integral on T e is zero because the integrand is analytic in the exterior of the unit 
disc and has a double zero at oo. The limit of the integral on / j e is given by 

1 f f t-w \ 2 &( T ) dT 1 



Therefore 



— lim 



lim 

e^o 27ri L \$(r)-$(w)/ r r - w 2w®'(w 



(61) — ^ = &(w)I 2 (w) + wGl 



= &(w)I 2 (w) + -\ 
dw w z 

Lemma [5] below applied for n = 2 yields I 2 G C^(T), for any < /3 < a and hence we 
conclude that w i— >■ belongs to the space C^(T), for any < /3 < a, that is, that 
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I\ G C 1+/3 (T). Thus from the expression (155|) . we find that 

(62) qeC 1+ P(T), 0e)O,a[. 

Before dealing with Lemma [5] we discuss how the smoothness of q translates into the 
same type of smoothness of More precisely, we will prove the following: for each 
positive integer n and < f3 < 1, 

(63) q G C n+/3 (T; T) =>- $ G C n+1+/3 (T; C). 

For this purpose, we first relate the regularity of the map w i— > q(w) to the smoothness of 
the Jordan curve $(T) and then we use the Kellogg- Warschawski Theorem |WSj to get 
a global regularity result for the conformal map $. 

Using the conformal parametrization 6 G R i— > $(e* 61 ), we get easily the following 
formula for the unit tangent vector f(9) to the curve $(T) at the point $(e* e ), 

Consequently 

(64) [t(6)} 2 = -w 2 ^), w = e ie . 

Since $ belongs to C 1+a } the map h-> f (0) must be in C a (R; T) C C(R; T) and by the 
lifting theorem there exists a continuous function : R — )■ R, such that 

f(0) = e im , 6 G R. 

Recall that we have established that q G C 1+/3 (T),0 < /3 < a and thus [r] 2 remains in 
the same space. Therefore, the argument function <p can be recovered by the formula 

HO) = HO) - / -Udt, with a(t) := [f(t)] 2 . 
2 Jo °W 

This shows that the argument function is in C 1+/3 (R) and consequently f G C 1+/3 (R). 
More generally, the preceding formula for shows that, for each non-negative integer n 
and p G]0,1[, 

q G C n+/3 (T; T)^fe C n+,3 (R; T). 
Now we will use the Kellogg- Warschawski theorem |WS] . which can be also found in jP| 
Theorem 3.6]. It asserts that if the boundary $(T) is a Jordan curve of class C n+1+I3 , 
with n a non-negative integer and < /3 < 1, then the conformal map $ : C\A — > C has 
a continuous extension to C\A which is of class C n+1+l3 . This means that 

t G C n+f5 {R; T) =^> $ G C n+1+/3 (T; C). 
Combining (162]) and fl63l) we obtain 

$ G C 2+/3 (T). 

We are now ready to iterate the preceding argument. Assume that $ G C" _1+/3 (T), 
< /3 < 1, for some n > 3. We are going to show that, 

(65) $GC n+7 (T), 7G]0,/3[. 

This will complete the proof that $ G C°°(T). 

We need the following general lemma. For $ G C n (T), let 

P B ($)(r,«;) = 53— f^(r-«;y 

3=0 J ' 

be the Taylor polynomial of degree n of $, around the point w, evaluated at the point r. 
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Lemma 5. Assume that $ G C n - 1+a (T), for n > 2, < a < 1. Let T n be the operator 
(66) 1 
with kernel 



T n g(w) = — / K n (w, r)g(r)dr, weT, g G L°°(T), 



K n (w,r 



$(r) - P n ^)(r,w) 



Then for any j3 g]0, a[ 



($(r) - $H) r 

^riflilc/^T) < C||5iU°°(T)- 



Proof. The lemma is easily proven by standard methods, as in |MOV} p. 419 ], once one 
knows that K„ satisfies 



(67) 
and 



\K n (w,T) < C\t- w\ a ~\ r,wGT,r^«; 



(68) \K n { Wl ,r)-K n {w 2 ,r)\<C 



\Wi - w 2 \ 



wi, u>2 G T, It — wi\ > 2\wi — w 2 \. 



\T — W\\ 

It is obvious that ( |67|) holds by Taylor's formula. For (168]) we write 



\K n {w U T) - K n (w 2 ,r)\ < 



($(r) - ^(w 1 )) n 



+ |$(r) - P n ^)(r,w 2 ) 
:= I + 11. 



($(r) - $(wi)) n ($(r) - ®(w 2 )Y 



The term II can easily be controlled via a gradient estimate by 

| II | < C\r-w 2 \ n - 
\wi - w 2 \ 



l+a 


Wi - W 2 






T — Wi 


n+1 



\T — Wi\ 

We have used in the last inequality the equivalence \\r — W\\ < \r — w 2 \ < §|t — W\\. The 
term / is estimated by observing that there is an elementary formula for the difference 
of two Taylor's polynomials around different points w± and w 2 , namely, 

n— 1 

(69) P n ^)(r,w 2 )-P n ^)(r, Wl ) = ^ ^(w 2 )-P n ^(^y)(w 2 , Wl ) 



(t - w 2 y 



3=0 



J- 



This follows easily from the identity, 



n-l 



P n ^)(T, Wl ) = J2 (T ^ {d^P n -l}(w 2 , Wl ) 



(r - w 2 ) j v-^ (w 2 - Wi) 



— iiu \ k 3 



-$( fc )l 



Formula fl69|) yields 



j=0 
n-l 

j=0 J ' 
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m < cy, 

j=0 



n-1 

\w 2 - W\\ 



l-j+a 


r - w 2 \ 


3 


1 






r — W\\ 


n 



< c 



\Wi - w 2 \ 
It — wA 



□ 



In view of fl53|) the only task left is to show that I\ G C n ~ 1+7 (T), < 7 < /3, provided 
$ G C n ~ 1+/3 (T). Indeed, this will lead to q G C n ~ 1+7 (T) and, according to the discussion 
above on the Kellog-Warschawski Theorem, we conclude that $ G C n+7 (T). Now, in 
order to prove that I x G C n ~ 1+7 (T) we need to establish a recursive formula for the 
higher order derivatives of Ii, which is the goal of the next lemma. 

Lemma 6. Let $ G C n+/3 (T), < < 1, n > 2, and set 



2m J T ($(r) - $(«;)) 



Then 



(70) = n$ , HWW,«ieT. 

aw 

Remark. Notice that formula (157]) for the derivative of I\ falls out of the scope of (1701) . 
Indeed, it is a fortunate fact that a compact formula as (T70|) can be found. 



The proof of the preceding lemma depends on the following calculation. 
Sublemma. Let $ G C X (T) andn>l. Then 

(71) p.*— / {T ~ WT +I nr)~ = t—\ ) WET. 

Proof. Consider again the closed Jordan curve T e and the arc j e used to deal with (150]) 
and (160]) . The principal value integral in the statement of the sublemma is the limit, as 
e tends to 0, of the sum of two terms. The first is the integral over T e of the integrand 
in ( 171 p . which is zero because the integrand is an analytic function of r in the exterior of 
the unit disc with a zero at 00 of order at least 2. The second term is minus the limit as 
e tends to of the integral of the same expression over the arc j e . Since 

1 f &(t) 1 1 
lim / -— , . dr = lim A~ log $ = - 

2m e^o 7 7e $( r ) - 2m e^o k 6 2 

is clear that the limit of the second term is — 1/(2 $>'(w) n w n ). Here A 7e log$ stands for 
the variation of log<3? on the arc 7 e . □ 



Proof of Lemma® Since $ G C n (T), the fraction in the integrand in /„ extends contin- 
uously to the diagonal of T taking the value (— l) n ^ n \w) w 2n /(n\ $'(u>) n ). We can then 
take derivatives inside the integral and the boundary terms arising in the integration by 
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parts are zero. Thus 



dw K ' 1 ,y 2m J T ($( r )_$( w )) rt+i 1 ^ 



+ 77 r^P-v. — / 7^ TT^n® ( T )dr 



n — l\ w 2 ' 2-ki J T ($(r) - $(u>)) 
T 1 (w) + T 2 (w). 



The next move consists in adding and subtracting $( n )(w)(r — w) n /n\ to the numerator 
of the fraction in the integrand of Ti(w). The result is 



^ H I \, P- v -— / 7T-^ . / ^n+i $ ( r ) rfr 



n-1! 2m y T ($( r ) _ 

$Wfi<7j 1 /■ f r - w) n 

= n ^ 1 W + ^)-iip.y.- / -J — pr ^(r)dr 

n— 1! ^« Jt ($(r) - $(»)) 

:= n^'H/n+iH+TsH 

We claim that T2(w) + T^w) = 0, which ends the proof of the lemma. The only 
difficulty is to compute the principal value integrals, which are the same except for a 
shift in the exponents of the integrand. For instance, by (|56|) and the sublemma we see 
that the principal value integral in the term T 3 (w) is 

1 f {T^) n {-l) n 1 f (r-w) n , dr 

p.v. / --r$ • (r)dr = - — p.v. / tt<& (t) — 

2m J T ($( r )-$H) n+1 w n 2m J T ($( r ) - $H) n+1 V 



1 



,n+l 



2[$'(w)] n w 2n ' 

This completes the proof of the Lemma □ 

We can now finish the proof of Theorem [21 Recall that our assumption is $ G 
C n - 1+ P(T) for some f3 g]0, 1[ and n > 3, and we want to conclude that h G C n " 1+7 (T), 
for any 7 G]0,/3[. We have already seen that this will give $ G C n+7 (T), 7 G]0,/3[. We 
apply Lemma [5] and we get that I n = T n ($') belongs to C 7 (T), < 7 < (3. Formula 
(1701) in Lemma |6j readily yields 7 n _i G C 1+7 (T). Iterating the use of (1701) we obtain 
h e C rt - 2+7 (T), and so, finally, by §2), h G C™- 1+7 (T). □ 

We end the paper with the following remark. 

Corollary. For each integer m > 3, there exists a small positive eo = €o{m) such that if 
f G V defines an m-fold V-state D, H/Hc^t) < e o an d \\f\\c 2 {r) < 1/2, then D is convex. 

Proof. As it is well-known, if D is a Jordan domain bounded by a smooth Jordan curve 
of class C 2 , then D is convex if and only if the curvature of the boundary curve does not 
change sign on the curve. By Theorem [2] we know that the boundary of our Estate D 
is of class C°°. To compute the curvature at boundary points we resort to the conformal 
parametrization $(e ie ). The velocity vector (or tangent vector) and the principal normal 
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to the curve at the point $>(w),w = e id , are given by 

v{6) = iw$ (w) and n{6) = 



\&(w)\ 

respectively. On the other hand, according to a well-known classical formula for the 
curvature k{9) at the point $(e* £l ), we have 



Ite(dg{$(e")}n(fl)) 
^1 = ^77n^ 



\v 



A straightforward computation yields 

d? g {$(e ie )} = -w(&(w)+w<S>"(w)). 

Thus the curvature is 

1 / $"W\ 

Since $(«;) = w + 

1 + W - 7 , = 1 + W f / > 

1 + /'(«;) 

and so 



i-l/'WI" i-|l/llc»m 

which is non-negative if ||/||c 2 (T) < 1/2. □ 

The reader will find interesting information on conformal mappings and convexity in 
Duren's book [Duj . 
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